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Abstract 

Measurements are reported of the production of dijet events with a leading neutron in ep 
interactions at HERA. Differential cross sections for photoproduction and deep inelastic 
scattering are presented as a function of several kinematic variables. Leading order QCD 
simulation programs are compared with the measurements. Models in which the real or 
virtual photon interacts with a parton of an exchanged pion are able to describe the data. 
Next-to-leading order perturbative QCD calculations based on pion exchange are found to 
be in good agreement with the measured cross sections. The fraction of leading neutron 
dijet events with respect to all dijet events is also determined. The dijet events with a 
leading neutron have a lower fraction of resolved photon processes than do the inclusive 
dijet data. 
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1 Introduction 



Previous HERA measurements [1,2] show that the cross section for the semi-inclusive ep scat- 
tering process 



where the leading neutron carries more than 70% of the proton beam energy, is reasonably well 
described by the pion exchange mechanism [3-7]. In this picture, the virtual photon interacts 
with a parton from the pion. Constraints on the pion structure function are thus obtained. How- 
ever, the Soft Colour Interaction model, in which colour neutral partonic subsystems are formed 
by non-perturbative soft gluon exchanges [8,9], describes the data equally well [1]. 

In the present analysis, the leading neutron production mechanism is investigated further by 
requiring that the system X in (1) contains two jets with large transverse momenta 



This allows more detailed comparisons of the measurements to be made with model predictions. 
In addition, the jet energy provides a hard scale which allows the comparison of perturbative 
QCD with the data for all photon virtualities Q 2 . The cross sections are measured in both pho- 
toproduction (Q 2 < 10~ 2 GeV 2 ) and deep inelastic scattering (DIS, 2 < Q 2 < 80 GeV 2 ). They 
are given as a function of Q 2 and of the kinematic variables of the jets. Monte Carlo predictions 
based on leading order (LO) QCD models are compared to the data, as are next-to-leading order 
(NLO) QCD calculations [10]. Furthermore, a detailed comparison of dijet production with and 
without the requirement of a leading neutron is made. In the photoproduction regime, similar 
studies have been reported by the ZEUS Collaboration [11]. 

2 Event Kinematics and Reconstruction 

The semi-inclusive reaction is sketched in Fig.Uk, in which the 4-vectors of the incoming 
and outgoing particles and of the exchanged photon are indicated. Figure [TJd depicts the dijet 
production reaction © under the assumption that it is mediated by pion exchange. 

The standard Lorentz invariant kinematic variables used to describe high energy ep interac- 
tions are the centre-of-mass energy squared s, the four-momentum transfer squared Q 2 and the 
inelasticity y: 



e + p — s-e + n + X, 



(1) 



e + p^-e + n + jet + jet + X. 



(2) 




s 



(k + P) 2 = 4:E e E } 



(3) 
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eft) 



P(P) 




Figure 1 : (a) A diagram for the process e+p — > e+n+X, (b) a diagram for the dijet production 
process e + p— > e + n + jet + jet + X assuming this proceeds via pion exchange. 



These are determined from the energies E e and E p of the lepton and proton beams, respectively, 
and from the energy E' e and polar angle 9' e of the scattered lepton in the laboratory frame 1 . 

Two more invariant variables, xl and t, are used to describe the kinematics of the semi-inclusive 
reaction <[T]): 

(q ■ p n ) E n 

x L — 



(q-P) E p 
2 ,~ _. 

XL ~' \ %L 



and t = (P -p n f ~ - (1 -a*) -mj) , (4) 



where P n is the neutron energy, p Tn is the momentum component of the neutron transverse 
to the direction of the incident proton and m n and m p are the neutron and proton masses, 
respectively. Experimentally, xl and t are determined from the measured energy and scattering 
angle of the leading neutron. 

In the pion exchange model, the photon interacts with a pion emitted from the proton. In 
this model, process (2) is represented by diagrams as sketched in Fig. [TJ5. The quantity x n 
denotes, neglecting masses, the fraction of the 4-momentum of the pion participating in the 
hard interaction 2 . It is related to x p , the fraction of the 4-momentum of the proton which enters 
the hard interaction, via x p = x n (l — xl). 

The quantity x 7 is the fraction of the 4-momentum of the photon which participates in the 
hard interaction. If the virtual photon is "resolved" and participates in the hard interaction via 
its partonic content, then s 7 < 1. If the interactions are "direct", i.e. the entire photon enters 
the hard scattering process, then x 7 = 1. 

'The right-handed HI coordinate system has its positive z direction along the proton beam direction and its 
origin at the nominal interaction point. 

2 The definitions of the variables Xl and are similar to the definitions of the variables (1 — xp) and (3, used 
in HI analyses of diffractive processes [12]. 
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The quantities x^ 1 , x 3 ^ 1 and x p et , which are estimators for x 7 , x w and x p , can be defined in 
dijet events using the jet transverse energies E^ and pseudorapidities r] jet according to: 



(E J T etl e-^ ea +E J T et2 e^ jet2 ) jet (E J T e fl e^ + E$*V eta ) 



2y.Be p 2E p 

and - g^!!±<i!!j (5) 

The pseudorapidity is defined by 77 = — In (tan |), where is the polar angle with respect to 
the z axis. 



3 Experimental Procedure 
3.1 HI detector 

The data used in this analysis were collected with the HI detector at HERA in the years 1996- 
97 and correspond to an integrated luminosity of 19.2 pb -1 . In these years the HERA collider 
was operated at positron and proton beam energies of E e = 27.6 GeV and E p = 820 GeV, 
respectively. 

A detailed description of the HI detector can be found elsewhere [13]. Here only the com- 
ponents relevant for the present measurement are briefly described. 

The e + p luminosity is determined with a precision of 1 .6% by detecting photons from the 
Bethe-Heitler process e + p — > e + p7 in the photon detector located at z — —103 m. The 
electron calorimeter of the luminosity system, located at z — —33 m, is used to detect the 
positrons scattered through very small angles (i.e. large 9' e ) in photoproduction events. Both 
detectors are TICl/TlBr crystal Cerenkov calorimeters with an energy resolution of a(E)/E = 
22%/y/E[GeV\. 

The Central Tracking Detector (CTD), with a polar angle coverage of 25° < 9 < 155°, 
is used to measure the trajectories of charged particles and to reconstruct the interaction ver- 
tex. The tracking system is surrounded by the finely segmented Liquid Argon (LAr) calorime- 
ter [14], which covers the range in polar angle 4° < 9 < 154° with full azimuthal accep- 
tance. The LAr calorimeter consists of an electromagnetic section with lead as absorber, and a 
hadronic section with steel as absorber. The total depth of the LAr calorimeter ranges from 4.5 
to 8 hadronic interaction lengths. Its energy resolution, determined in test beam measurements, 
is cr(E)/E w 12%/y^GeV] © 1% for electrons and a{E)/E w 50%/y/E [GeV] © 2% for 
charged pions. The absolute electromagnetic energy scale is known with a precision of 1%. 
The absolute hadronic energy scale for the jet energies used in this analysis is known with a 
precision of 4%. 

The polar angle region 153° < 9 < 177.8° is covered by a lead/scintillating-fibre calorime- 
ter, the SpaCal [15], with both electromagnetic and hadronic sections. For positrons, the energy 
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resolution is a(E)/E « 7.1%/ ^/ E[GeV] © 1% and the energy scale uncertainty is less than 
2%. The positron polar angle is measured with a precision of 1 mrad. The hadronic energy 
scale in the SpaCal is known with a precision of 7%. A Backward Drift Chamber (BDC) in 
front of the SpaCal is used to suppress background from neutral particles which can fake the 
scattered positron signal. 

The LAr and SpaCal calorimeters are surrounded by a superconducting solenoid which 
provides a uniform magnetic field of 1.15 T along the beam direction. 

Leading neutrons are detected in the Forward Neutron Calorimeter (FNC), which consists 
of interleaved layers of 2 m long lead strips and scintillator fibres and is located 107 m away 
from the nominal HI interaction point in the proton beam direction (for details see [1]). The 
energy resolution of the calorimeter is a(E)/E ps 20% for neutron energies between 300 and 
820 GeV. The absolute hadronic energy scale is known with a precision of 5%. Two segmented 
planes of hodoscopes situated in front of the FNC are used to veto charged particles. Each 
plane is constructed of 1 cm thick hexagonal scintillator tiles, which have the same lateral size 
as the calorimeter modules. The neutron detection efficiency of the FNC is (93 ± 5)%, the 
losses being due to the back- scattering of charged particles from the hadronic shower caused by 
the neutron which give signals in the veto hodoscopes. The acceptance of the FNC is defined 
by the aperture of the HERA beam line magnets and is limited to neutron scattering angles of 
9 n < 0.8 mrad with approximately 30% azimuthal coverage. 

3.2 Event selection 

The events used in this analysis are triggered by the coincidence of a track in the CTD with 
an electromagnetic cluster either in the SpaCal (DIS sample) or in the electron calorimeter of 
the luminosity system (photoproduction sample). A number of selection criteria are applied in 
order to suppress background and to confine the measurements to those regions of phase space 
where the acceptance is large and uniform. 

The reconstructed z coordinate of the event vertex is required to be within ±30 cm (~ 3a) 
of the mean z position of the interaction point. 

In the photoproduction analysis, the scattered positron energy E' e is limited to the range 
0.3 < y « (1—E' e /E e ) < 0.65. This cut is defined by the geometrical acceptance of the electron 
calorimeter which also restricts the value of Q 2 to be less than 10~ 2 GeV 2 . To ensure that the 
effect of radiative corrections in photoproduction is small and to suppress events in which a 
Bethe-Heitler event and a photoproduction event are superimposed, the energy measured in the 
photon detector of the luminosity system is required to be less than 1.5 GeV. 

The final state positron candidates in DIS are required to have polar scattering angles in the 
range 156° < 9' e < 176° and energies E' e greater than 10 GeV. The energy and angle, E' e and 
9' e , are determined from the associated SpaCal cluster in combination with the interaction vertex 
reconstructed in the Central Tracker. The analysis is restricted to the region 2 < Q 2 < 80 GeV 2 
and 0.1 < y < 0.7. To suppress events with initial state hard photon radiation, as well as 
events originating from non-ep interactions, the quantity E—p z , summed over all reconstructed 
particles including the positron, is required to lie between 35 GeV and 75 GeV. This quantity, 
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which refers to the energy and longitudinal momentum component of each final state particle, 
is expected to be twice the electron beam energy for contained events. 

All events that satisfy the selection cuts are subjected to a jet search using a cone algorithm 
with radius R = 1 [16]. The jet finding is performed using the energies measured in the LAr 
and SpaCal calorimeters in the 7*p centre-of-mass frame, with transverse energies calculated 
relative to the 7* axis in that frame. To ensure that the bulk of the jet energy is contained within 
the LAr calorimeter, the laboratory pseudorapidity of each jet axis is restricted to the region 
— 1 < rjl^l < 2. Events which have exactly two jets are selected. The transverse energies of 
these jets must be above 7 GeV (first jet) and above 6 GeV (second jet). 

For the cross section measurements, events with a leading neutron are selected from the 
inclusive dijet samples by requiring a cluster in the FNC with an energy above 500 GeV. At 
such high energies the background contribution due to other neutral particles is negligible. 

The final photoproduction data sample contains about 69 000 dijet events, of which 372 
events contain a leading neutron with energy E FNC > 500 GeV. In the DIS sample 23 000 
dijet events are selected, of which 213 satisfy the neutron identification criteria with -Efnc > 
500 GeV. 

The kinematic regions within which this measurement is made are summarized in Table 1. 





Kinematic regions 


Photoproduction 
DIS 


Q 2 < 10~ 2 GeV 2 , 0.3 < y < 0.65 
2 < Q 2 < 80 GeV 2 , O.K y < 0.7 


Dijets 
Neutrons 


E 3 T etl > 7 GeV, E 3 T et2 > 6 GeV, -1 < rf^' 2 < 2 
E n > 500 GeV, 9 n < 0.8 mrad 



Table 1 : The kinematic regions within which the cross sections are measured. 



4 Monte Carlo Models 

Monte Carlo samples are used to correct the data for inefficiencies, acceptance effects, migra- 
tions and the effects of QED radiation. They are also used to correct for hadronization effects in 
the comparison with NLO QCD calculations. Monte Carlo predictions based on several leading 
order QCD models are also compared with the data. 

All the Monte Carlo programs generate hard parton-level interactions using the Born level 
QCD matrix elements with a minimum cut-off on the transverse momentum of the outgoing 
partons. They differ in the assumptions made about the origin of the partons, which may come 
from either the incident proton or an exchanged pion, and in the details of the hadronization 
models. After hadronization, the response of the HI detector to the events is simulated in detail 
and they are passed through the same analysis chain as is used for the data. 

In addition to the models which are frequently used in studies of inclusive jet production, 
namely PYTHIA [17] for photoproduction and RAPGAP [18] and LEPTO [8] for DIS, models 
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in which the hard interaction proceeds only via 7r-exchange are also used: the 7r-exchange ver- 
sion of RAPGAP for both photoproduction and DIS and POMPYT [19] for photoproduction. A 
model in which a colour neutral system is formed non-perturbatively by soft colour interactions 
(SCI) [9] is also compared with the data. This mechanism is implemented in the Monte Carlo 
program LEPTO. 

The PYTHIA event generator simulates hard photon-proton interactions via resolved and 
direct photon processes. It is used with a minimum value for the transverse momenta of the 
outgoing partons in the hard interaction process (j)™ m ) of 2 GeV. The GRV-LO parton den- 
sities are used for the photon [20] and the proton [21]. The photon flux is calculated in the 
Weizsacker- Williams approximation [22]. Higher order QCD radiation effects are simulated 
using initial and final state parton showers in the leading log approximation. The subsequent 
fragmentation follows the Lund string model as implemented in JETSET 7.4 [23]. PYTHIA 
can also simulate multi-parton interactions (MI), which are calculated as LO QCD processes 
between partons from the remnants of the proton and the resolved photon. The resulting ad- 
ditional final state partons are required to have transverse momenta above a cut-off value of 
1.2 GeV. It has previously been shown [24] that these additional interactions improve consid- 
erably the description of inclusive jet photoproduction. This option of PYTHIA is referred to as 
PYTHIA-MI below. The PYTHIA calculation is performed with version 5.7 and cross-checked 
with version 6.1. 

The program LEPTO 6.5 [8] generates DIS events. It is based on leading order electroweak 
cross sections and takes QCD effects into account to order a s . As in PYTHIA, higher order 
QCD effects are simulated using leading log parton showers and the final state hadrons are 
obtained via Lund string fragmentation. Higher order electroweak processes are simulated us- 
ing DJANGO [25], an interface between LEPTO and HERACLES [26]. The LEPTO program 
allows the simulation of soft colour interactions [9], through which the production of leading 
baryons and diffraction-like configurations is enhanced via non-perturbative colour rearrange- 
ments between the outgoing partons. In the following, the predictions based on this approach 
are denoted LEPTO-SCI. 

The program RAPGAP 2.8 [18] is a general purpose event generator for inclusive and 
diffractive ep interactions. In DIS, the RAPGAP simulation includes a contribution from re- 
solved photon events in which the photon structure is parameterized according to the SaS- 
2D [27] parton densities. These give a reasonable description of inclusive dijet production at 
low Q 2 [28]. In the version denoted below as RAPGAP- ir, the program simulates exclusively 
the scattering of virtual or real photons off an exchanged pion. Here, the cross section for 
photon-proton scattering to the final state nX takes the form 

where f n + / p (x L , i) is the pion flux associated with the beam proton and rf<T 7 * 7r+ ^ x denotes the 
cross section for the hard photon-pion interaction. The pion flux factor is taken from [4]: 
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where m n is the pion mass, g^/Air = 13.6 is the pirn coupling constant, known from phe- 
nomenological analyses of low-energy data [29], and R wn = 0.93 GeV -1 is the radius of the 
pion-neutron Fock state of the proton [4]. For the range of t and x L relevant here, this flux 
parameterization yields results very similar to those of the parameterization of [5], which in- 
cludes the full t dependence expected from the pion Regge trajectory and which is used in 
our previous analysis [1]. If not otherwise stated, the GRV-7T-LO [30] pion structure function 
parameterization is used. 

The POMPYT 2.6 [19] Monte Carlo program is an extension of PYTHIA, which models 
colour singlet exchange processes in photoproduction. For pion exchange processes, POMPYT 
simulates the scattering of real photons off the exchanged pion, using the pion flux parameteri- 
zation of eq.(7). POMPYT yields results very similar to those of RAPGAP-7T. 



5 Systematic Uncertainties 

The acceptance of the FNC calorimeter is defined by the interaction point and the geometry 
of the beam guiding magnets and is determined using Monte Carlo simulations. The angular 
distribution of the neutrons produced in reaction (1) is sharply peaked in the forward direction, 
and the observed cross section therefore depends critically on small inclinations of the incoming 
proton beam with respect to its nominal direction. This effect is studied with the help of the 
Monte Carlo simulations described above. The overall uncertainty in the FNC acceptance is 
estimated to be about 10%. The uncertainty in the neutron detection efficiency leads to an 
additional 5% systematic error, and the uncertainty of the FNC absolute energy scale to a 6% 
systematic error. These effects contribute to the overall normalization error. 

The 4% uncertainty on the absolute hadronic energy scale of the LAr calorimeter leads to 
an uncertainty of about 15% on the jet cross section. This is strongly correlated between data 
points. 

The uncertainty of the acceptance of the electron calorimeter of the luminosity system, in 
the photoproduction case, amounts to about 6%. This includes the uncertainty on the luminosity 
measurement of 1 .6% and contributes to the overall normalization error. 

The uncertainties on the measurements of the positron energy and angle in the SpaCal lead 
to 6% systematic uncertainties in the DIS cross sections. 

As shown below, the models based on the pion-exchange mechanism describe the data well 
and are therefore used to estimate acceptance and migration corrections. These corrections are 
determined from the POMPYT Monte Carlo simulations in the photoproduction case and using 
RAPGAP-7T for DIS. The uncertainties of the corrections are estimated from the differences in 
the results when other models are used: PYTHIA or RAPGAP-7T in the photoproduction case 
and LEPTO or RAPGAP in the DIS case. The estimated uncertainties are between 10% and 
15% for all distributions. 

Due to the energy cut in the photon detector, QED radiative corrections are small [31] in the 
photoproduction case and are neglected here. For the DIS sample, the QED radiative corrections 
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amount to less than 10%, as evaluated using RAPGAP interfaced to the HERACLES program. 
The uncertainty arising from the radiative corrections is about 5%. 

Finally, a 3% correlated uncertainty is attributed to the trigger efficiencies as evaluated using 
other, independent triggers. 

The correlated error contributions are only weakly dependent on the kinematic variables 
studied, causing a normalization uncertainty of about 20% on the cross sections in both the 
DIS and photoproduction cases. The uncorrected, point-to-point systematic uncertainties range 
from 11% to 17%. 

In the figures, the outer error bars represent the quadratic sum of the point-to-point sys- 
tematic errors and the statistical errors, while the inner error bars show the statistical errors. 
The normalization error is not shown in the figures, but is included in the tables as a correlated 
systematic error. 

6 Results and Discussion 

6.1 Neutron energy distribution 

The neutron energy spectrum allows discrimination between different models for the production 
of leading neutron events. Figure |2] shows the energy spectra for the photoproduction and DIS 
dijet samples, as measured in the FNC, for energies above 400 GeV. The data are not corrected 
for efficiencies, acceptance or migration between bins. These effects, however, are taken into 
account in the Monte Carlo simulations. In Fig. Eh the photoproduction and DIS data are shown, 
normalized to the total number of events in the corresponding inclusive dijet samples. There is 
a significantly higher fraction of leading neutrons in DIS dijet events than in photoproduction 
dijet events. However, the shape of the energy spectrum is similar in both samples. 

In Figs. Eb-d the photoproduction and DIS data are shown together with the predictions 
from several Monte Carlo models. The model predictions are all normalized to the integrated 
luminosity of the corresponding data samples. 

For simulated neutron energies E n > 400 GeV, the photoproduction data in Fig.Eb are rea- 
sonably well described in shape and magnitude by the 7r-exchange model RAPGAP- it, as well 
as by PYTHIA without multiple interactions. If multiple interactions are included, PYTHIA 
fails to describe the data, predicting a rate which is too high for neutron energies between 400 
and 600 GeV. This observation is in contrast to the inclusive jet measurement (without the re- 
quirement of a leading neutron), which is described by PYTHIA only if multiple interactions 
are included, especially at low jet E T [24]. 

The RAPGAP predictions are compared with DIS data in Fig.|2j;. The 7r-exchange version 
describes the shape of the distribution well, but somewhat overestimates the absolute rate, while 
the rate predicted by the standard RAPGAP DIS version is too low. The LEPTO prediction 
is also too low, as is shown in Fig. EJi. The LEPTO predictions are somewhat increased for 
E n £ 500 GeV if soft colour interactions (LEPTO-SCI) are included. 

To summarize, the pion exchange models describe the shape of the observed FNC energy 
spectra well for energies above 500 GeV, as does LEPTO with soft colour interactions. 
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Figure 2: Energy distributions observed in the FNC for dijet events: (a) Photoproduction data 
compared with DIS data, normalized to the respective number of dijet events in each inclusive 
dijet sample; (b) The Monte Carlo model predictions for RAPGAP- ir, PYTHIA and PYTHIA- 
MI, compared with photoproduction data; (c, d) The Monte Carlo models RAPGAP (both the 
7r-exchange and the standard DIS versions) and LEPTO (with and without SCI), compared with 
DIS data. The errors on the data points are statistical only. The corresponding kinematic regions 
for the photoproduction, DIS and dijet selections are given in Table 1 . 
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6.2 Cross section measurements 

In this section, differential cross sections at the hadron level are presented for dijet production in 
the photoproduction and DIS regimes for neutron energies E n > 500 GeV, corresponding to the 
region in which 7r-exchange models give a good description of the -Efnc distribution. The data 
are corrected for detector inefficiencies and migrations due to detector resolution effects using 
the Monte Carlo simulations described in sections 4 and 5. The results are given in Figures 3-6 
and Tables 2-3. 

In Fig.|3]the jet cross sections are shown as a function of EjT and rf^ h for the photoproduc- 
tion and DIS regimes using both jets in the event. 

Taking the 20% normalization uncertainty into account, the data are well described by the 
7r-exchange model RAPGAP-7T in both DIS and photoproduction. However, PYTHIA without 
multiple interactions, which does not include pion exchange, also provides a good description 
of the photoproduction data. The inclusion of multiple interactions in PYTHIA results in a 
predicted cross section which is too high for low values of E^ 1 and for values of rf^ b > —0.5. 
It is also seen from Figs. OJ? and|3ji that standard DIS processes, as simulated by the LEPTO 
program, tend to lie below the data, even if soft colour interactions are included. 

The measured Q 2 dependence of the dijet DIS cross section is shown in Fig.|3 together with 
predictions from Monte Carlo simulations. The 7r-exchange version of RAPGAP describes 
the measured distribution fairly well, whereas LEPTO with SCI reproduces the shape of the 
distribution but yields a lower prediction over the whole Q 2 range. 

The dependence of the dijet cross section on the fractional momenta x^ 1 and x{ et , deter- 
mined according to eq.(5), is shown in Figs. |5] and The measured xif 1 distribution in the 
photoproduction regime, together with the RAPGAP- 7r and PYTHIA model predictions shown 
in Fig.|5k, clearly demonstrate the large contribution of resolved photon processes. The shape 
of the distribution is well described by PYTHIA and RAPGAP-7T. PYTHIA with multiple in- 
teractions predicts too high a cross section at x^ 1 < 0.6. In the DIS regime, as is clear from 
the x^ 1 distribution shown in Fig. |5J), direct photon interactions dominate. However, a small 
fraction (~ 15%) of resolved photon interactions is necessary to fully describe the data with 
the RAPGAP-7T simulation. The LEPTO-SCI model, which does not include resolved photon 
processes, provides a poor description of the shape of the distributions. 

The x{ e * distributions in the photoproduction and the DIS samples, shown in Figs. [6^ and 
03, respectively, are similar in shape. The 7r-exchange model RAPGAP- n, and the PYTHIA 
and LEPTO-SCI models, provide a fair description of the data. 

For the RAPGAP- n predictions shown in Figs. 3 to 6 the GRV-7T-LO parameterization of 
the pion parton distribution functions (PDFs) is used. However, within the 20% normalization 
uncertainty, a similar quality of description is provided if other parameterizations of the pion 
PDFs [32,33] are used. 

6.3 Comparison of the photoproduction cross section with NLO QCD 

The parton level cross sections for dijet photoproduction in photon-pion and photon-proton 
collisions are calculated in QCD at both leading and next-to-leading order [10]. In the latter 
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Figure 3: The measured differential ep cross sections as a function of E 3 ^ 1 and "i]\f h for dijet 
events with a leading neutron. The cross section is given for photoproduction (a,c) and DIS (b,d) 
and compared with the Monte Carlo simulations. Inner error bars show the statistical errors, 
while the outer error bars represent the statistical and systematic errors, added in quadrature. 
The overall normalization uncertainty of 20% is not shown. The kinematic regions within which 
this measurement is made are given in Table 1 . 



case, infrared and collinear singularities are cancelled using the phase space slicing method 
with an invariant mass cut-off. The renormalization and factorization scales are defined to be 
the maximum transverse energy of the outgoing partons. The scale uncertainty amounts to 
approximately 15% on average, but is significantly larger (up to 30%) for low x J y et and high 
rjjal, as estimated by varying the scales by factors of 0.5 and 2. The photon flux is calculated 
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Figure 4: The measured deep inelastic ep cross section as a function of Q 2 for dijet events 
with a leading neutron. Inner error bars show the statistical errors, while the outer error bars 
represent the statistical and systematic errors, added in quadrature. The overall normalization 
uncertainty of 20% is not shown. Predictions from Monte Carlo simulations are compared with 
the measurements. The kinematic regions within which this measurement is made are given in 
Table 1. 

using the Weizsacker- Williams approximation [22]. The light-cone form factor [4] is used in 
the pion flux with the same parameters as for the RAPGAP-7T Monte Carlo predictions. In the 
calculations, the GRV parameterizations for the parton distribution functions are used for both 
the photon and the pion. A cone algorithm with radius R = 1 is used in the definition of jets. 

Since the QCD calculations refer to jets of partons, whereas the measurements refer to jets 
of hadrons, the predicted cross sections are corrected to the hadron level using factors evaluated 
from the LO Monte Carlo programs described in section 4. The hadronization correction factor, 
(1 + Shadr), is defined as the ratio of the cross section obtained with jets reconstructed from 
hadrons to that using jets reconstructed at the parton level after the generation of parton showers. 
The corrections are calculated by taking an average of the results from two different Monte 
Carlo models (POMPYT and RAPGAP-7r). The uncertainty of these corrections is taken to 
be half the difference between the results obtained from the two models, which is typically 
smaller than 5%. The hadronization corrections have a tendency to increase the calculated 
NLO cross section at low (by approximately 30% for E^ 1 = 7 GeV) and to decrease the 
cross section at high E^ 1 (by approximately —10% for E^ = 20 GeV). As a function of 
xtf 1 , the hadronization corrections increase the cross section by about 25% at the lowest xt^ 1 , 
by about 100% for x^ 1 = 0.8, and are close to zero for x^ 1 = 1. The hadronization corrections 
show only a weak dependence on rf^ and x{ e *. 

The LO and NLO calculations are shown with the measured cross sections in Fig. [71 After 
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Figure 5: The measured differential ep cross sections as a function of x^ 1 for dijet events 
with a leading neutron. The cross section is given for photoproduction (a) and DIS (b) and 
compared with the Monte Carlo simulations. Inner error bars show the statistical errors, while 
the outer error bars represent the statistical and systematic errors, added in quadrature. The 
overall normalization uncertainty of 20% is not shown. The kinematic regions within which 
this measurement is made are given in Table 1 . 



corrections for hadronization, there is good agreement between the NLO calculations and the 
measurements. The LO and NLO predictions without hadronization corrections do not describe 
the data. 



6.4 Ratios of leading neutron to inclusive dijet cross sections 

The ratio of the dijet cross sections with and without the requirement of a leading neutron, f in , 
is an interesting quantity in that it discriminates between the various Monte Carlo models used 
to describe leading neutron production. Further, if the hard interaction is independent of the 
neutron production, it should be essentially independent of the jet kinematics which reflect the 
hard process, neglecting possible phase space effects. As the only difference in the event selec- 
tion for the leading neutron data and the inclusive dijet samples is the requirement of a leading 
neutron, some important systematic uncertainties cancel in the ratio. The remaining overall 
normalization uncertainties, mainly associated with the acceptance, efficiency and energy scale 
of the FNC calorimeter, are about 13%. 

For the photoproduction data, //„ is shown in Fig. [5] and in Table 4 as a function of the jet 
variables. Figure shows that j\ n is, within errors, independent of and has an average 
value of about 2.3%. However, as can be seen in Figs.Hb-d, fi n shows a dependence on r\\f h , 
x J 7 et and x J p et . These dependences can only partly be reproduced by the PYTHIA model, which 
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Figure 6: The measured differential ep cross sections as a function of x^ 1 for dijet events 
with a leading neutron. The cross section is given for photoproduction (a) and DIS (b) and 
compared with the Monte Carlo simulations. Inner error bars show the statistical errors, while 
the outer error bars represent the statistical and systematic errors, added in quadrature. The 
overall normalization uncertainty of 20% is not shown. The kinematic regions within which 
this measurement is made are given in Table 1 . 

provides some estimate of the size of possible phase space effects. A better description of the 
ratio in Fig. [5] is possible, if the leading neutron data are described by the 7r-exchange model, 
RAPGAP-7T, and the inclusive dijet data by PYTHIA-MI. This comparison suggests that the 
mechanism for dijet production in events with a leading neutron differs from that in inclusive 
dijet events. 

This is further studied by measuring the Q 2 dependence of the ratio of the DIS dijet cross 
sections with and without the leading neutron requirement. The result is shown in Fig. |9] and 
in Table 5. Here, the point at Q 2 = is the average of the ratios for photoproduction shown in 
Fig.[5J Within the experimental errors, the RAPGAP model describes the measured ratio, when 
the leading neutron and the inclusive dijet data are represented by RAPGAP- n and standard 
RAPGAP, respectively. However, there is some tendency for the measured ratio to increase with 
Q 2 , for Q 2 below 20 GeV 2 . A similar Q 2 dependence was observed by the ZEUS Collaboration 
in the analysis of inclusive DIS events with leading neutrons [2]. 



6.5 Discussion 

It was observed in [1] that pion exchange provides a good description of the semi-inclusive 
DIS process (1) in which a leading neutron is produced. The present results, given in sections 
6.1 to 6.3, demonstrate that this is also the case for the small subsample of leading neutron 
events in which a dijet system is produced, in both DIS and photoproduction. This observation 
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Figure 7: The measured differential ep photoproduction cross section as a function of , 
Vial* an d x i et f° r dijet events with a leading neutron. Inner error bars show the statisti- 
cal errors, while the outer error bars represent the statistical and systematic errors, added in 
quadrature. The overall normalization uncertainty of 20% is not shown. LO, and NLO QCD 
predictions [10] before and after correction for hadronization effects, are compared with the 
measurements. The shaded bands show the quadratic sum of the factorization and renormaliza- 
tion scale uncertainties of the NLO predictions and of the uncertainty due to the hadronization 
corrections, 1 + Shadr- The kinematic regions within which this measurement is made are given 
in Table 1 . 
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Figure 8: The ratio of the cross section for dijet photoproduction with a leading neutron to 
that for inclusive dijet photoproduction, as a function of E 3 ^ 1 , r] 3 ^, x^ 1 and x 3 p et . Inner error 
bars show the statistical errors, while the outer error bars represent the statistical and systematic 
errors, added in quadrature. The overall normalization uncertainty of 13% is not shown. Monte 
Carlo predictions for the ratios are obtained by using either RAPGAP-7T for the leading neutron 
cross sections and PYTHIA-MI for the inclusive cross sections, or by using PYTHIA in both 
cases. The PYTHIA prediction without multiple interactions has been scaled by a factor 0.6 to 
ease the shape comparison. The kinematic regions within which this measurement is made are 
given in Table 1 . 
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the outer error bars represent the statistical and systematic errors, added in quadrature. The 
overall normalization uncertainty of 13% is not shown. The Monte Carlo prediction, shown only 
for DIS, is obtained by using RAPGAP-7T for the leading neutron cross section and RAPGAP 
for the inclusive cross section. The kinematic regions within which this measurement is made 
are given in Table 1. 

is not trivial, as the parameters used in the Monte Carlo models to empirically describe the 
pion exchange were determined in hadronic reactions and no tuning to the present data was 
performed. 

It is also observed that, for the DIS sample, the standard Monte Carlo models for the sim- 
ulation of the hadronic final state, such as LEPTO and RAPGAP, predict cross sections for the 
production of dijets with a leading neutron which are too small. Again, no attempt to tune pa- 
rameters was made. The increase of the leading neutron rate in dijet production caused by the 
introduction of non-perturbative soft colour interactions in LEPTO-SCI is not large enough to 
provide a good description of the measurements. However, LEPTO-SCI successfully describes 
the DIS reaction ep — > enX [1]. 

For the photoproduction sample, equally good descriptions are obtained with the pion ex- 
change model RAPGAP-7T and the standard Monte Carlo program PYTHIA. The predictions 
of PYTHIA-MI clearly fail to describe the leading neutron data. However, the introduction of 
multiple interactions in PYTHIA is necessary to describe inclusive jet production [24]. The x^ 1 
distributions in Fig. demonstrate that this discrepancy between the data and PYTHIA-MI 
predictions is due to the poor description of resolved photon processes. The relative fraction of 
these processes is considerably lower in dijet production with a leading neutron than in inclu- 
sive dijet production. This finding is corroborated by the ratios of the dijet cross sections with 
and without the leading neutron requirement, presented in Fig. Ht as a function of x^ 1 . This 
ratio increases by a factor of 3 as x^ et increases from 0.3 to 0.9. The tendency seen in Fig. 



CM 

-„ 0.045 

SP 0.04 
D 

II 0.035 
a 

0.03 
0.025 

0.02 
0.015 

0.01 
0.005 



20 



for this ratio to rise with Q 2 may be due to absorptive corrections, as has been pointed out by 
several authors [34] . 

The differences between the dijet production data with and without leading neutrons, as 
well as the kinematic dependences of the cross section ratios, shown in Figs. [HJ5— d, point to 
differences in the production mechanism of events with and without leading neutrons. The 
present analysis shows that pion exchange is able to describe the properties of leading neutron 
events. 

A similar analysis in the photoproduction regime has been published by the ZEUS Col- 
laboration [11]. Qualitatively, there is good agreement between the two studies. A detailed 
comparison is difficult, however, since the kinematic ranges and the jet algorithm used are dif- 
ferent. The ZEUS Collaboration also observes that the ratio of leading neutron events to the 
inclusive sample is independent of £ff*, and sees a similar dependence of this ratio on x^ 1 . The 
dependence observed in [1 1] as a function of r/ je * is somewhat weaker than that in Fig.[SJx 



7 Summary 

The production of dijet events with a leading neutron of energy E n > 500 GeV and polar 
angle 6 n < 0.8 mrad is studied in photoproduction (Q 2 < 10~ 2 GeV 2 and 0.3 < y < 0.65) 
and in deep inelastic scattering (2 < Q 2 < 80 GeV 2 and 0.1 < y < 0.7). Dijet events with 
E J T etl > 7 GeV and E J T et2 > 6 GeV are selected using a cone algorithm in the 7*p frame. The 
laboratory pseudorapidities of the jets are restricted to the region — 1 < rf^' 2 < 2. Differential 
cross sections are presented as a function of E^ 1 , rf^, x^ 1 and x{ e * for photoproduction, and as 
a function of Q 2 , E^*, rf^jj, x^ 1 and xl et for DIS. 

Both the cross section measurements and the neutron energy spectrum are reasonably well 
described by pion exchange models in which the photon interacts with a parton from the ex- 
changed pion. These models are based on leading order QCD. The phenomenological param- 
eters describing the pion exchange are taken from previous analyses of hadronic interactions. 
Next-to-leading order QCD calculations, after corrections for hadronization effects, also de- 
scribe the measured photoproduction dijet distributions, in normalization as well as in shape. 
The experimental uncertainties are still too large to discriminate in this kinematic region be- 
tween different parameterizations of the pion parton densities. 

Monte Carlo programs which are not based on the pion exchange mechanism, such as 
PYTHIA, are also able to describe the leading neutron dijet photoproduction data. However, 
the predictions of PYTHIA including multiple interactions fail to describe the measurements. 
For the DIS sample, the standard LO Monte Carlo models RAPGAP and LEPTO, which do 
not include meson exchange, give a poor description of the data. This remains true when soft 
colour interactions are added to LEPTO, in which case, however, the inclusive leading neutron 
data are described. 

The ratios of the cross sections with and without the leading neutron requirement are stud- 
ied as a function of Q 2 and of the jet kinematic variables listed above. If the hard process is 
independent of the leading neutron production, these ratios should not depend either on Q 2 or 
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the jet variables. Indeed, there is no evidence for a strong dependence of the ratio on Q 2 , or on 
E^* in photoproduction. However, the ratio in photoproduction rises with x^ 1 . This suggests 
that the leading neutron dijet data have a lower fraction of resolved photon processes than do 
the inclusive dijet data. 
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Table 2: The differential ep photoproduction cross section as a function of E^, rf^, xi^ and 
x jet f or ^ygj even j s Yyjjh a leading neutron. The kinematic regions within which this measure- 
ment is made are given in Table 1 . 
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Table 3: The differential deep inelastic ep cross section as a function of rjj^, xif 1 , x{ et and 
Q 2 for dijet events with a leading neutron. The kinematic regions within which this measure- 
ment is made are given in Table 1. 
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Table 4: The ratio of the cross section for dijet photoproduction with a leading neutron to that 

The kinematic 



for inclusive dijet photoproduction, as a function of E; 
regions within which this measurement is made are given in Table 1. 



•jet 
T 



jet 



x J y et and x J p et . 



Q 2 

[GeV 2 ] 


fin 


0~stat. 


Ouncorrel.syst. 


ficorrel.syst. 




< io- 2 




0.0233 


0.0012 


0.0016 


0.0029 


2 




6 


0.0241 


0.0029 


0.0017 


0.0031 


6 




12 


0.0282 


0.0041 


0.0020 


0.0061 


12 




20 


0.0320 


0.0053 


0.0022 


0.0041 


20 




40 


0.0307 


0.0050 


0.0021 


0.0039 


40 




80 


0.0279 


0.0056 


0.0020 


0.0036 



Table 5: The ratio of the cross section for dijet production with a leading neutron to that for 
inclusive dijet production, as a function of Q 2 . The kinematic regions within which this mea- 
surement is made are given in Table 1 . 
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